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We propose a detailed investigation of the schematic mode-coupling approach to experimental 
data, a method based on the use of simple mode-coupling equations to analyze the dynamics of 
supercooled liquids. Our aim here is to clarify different aspects of this approach that appeared so 
far uncontrolled or arbitrary, and to validate the results obtained from previous works. Analyzing 
the theoretical foundations of the approach, we first identify the parameters of the theory playing 
a key role and obtain simple requirements to be met by a schematic model for its use in this 
context. Then we compare the results obtained from the schematic analysis of a given set of 
experimental data with a variety of models and show that they are all perfectly consistent. A 
number of potential biases in the method are identified and ruled out by the choice of appropriate 
models. Finally, reference spectra computed from the mode-coupling theory for a model simple 
liquid are analyzed along the same lines as experimental data, allowing us to show that, despite the 
strong simplification in the description of the dynamics it involves, the method is free from spurious 
artifacts and provides accurate estimates of important parameters of the theory. The only exception 
is the exponent parameter, the evaluation of which is hindered, as for other methods, by corrections 
to the asymptotic laws of the theory present when the dynamics is known only in a limited time or 
frequency range. 



I. INTRODUCTION 

Since the mode-coupling theory (MCT) of the liquid- 
glass transition has been proposed in the mid-eightiescra, 
considerable work has been done in order to compare its 
predictions with experimental data and computer simu- 
lation resultsEl. 

Most of these tests of the MCT are based on the use 
as fitting functions of the quasi-universal laws that are 
obtained from an asymptotic analysis of the family of 
non-linear equations to which the mode-coupling equa- 
tions belong and that are valid only in the vicinity of the 
predicted so-called ideal glass transition singularity. For 
experimental data, this is essentially the only practica- 
ble approach, since typical glass-formers are molecular 
liquids for which a first-principle study is out of reach 
and since moat experiments, like dynamic light-scattering 
spectroscopy!], probe complex mechanisms in which var- 
ious types of dynamical processes are entangled. In the 
case of computer simulations, the situation appears more 
comfortable, since simple systems can be considered for 
which all the data needed for an extensive comparison 
with the theory can be directly obtained with high pre- 
cision. However, it turns out that one still has to rely on 
asymptotic results of the theory at some point. Indeed, 
one finds that the MCT is not able to predict accurately 
the location of the putative ideal glass transition: it has 
thus to be determined from the simulation results, usu- 
ally with a fit to an asymptotic law, in order to-calibrate 
the comparison between theory and simulation^. 

Such a use of the asymptotic laws does not go with- 
out difficulties, originating in the fact that the ideal glass 
transition is never observed as such: it is usually argued 
to be avoided because of the existence of so-called "ac- 
tivated processes" not taken into account by the theory, 
whose effect is to restore the ergodicity of the system in 



the predicted ideal glassy state. Thus the transition is ex- 
pected to be replaced by a smooth crossover between two 
dynamical regimes, one dominated by the mode-coupling 
mechanism and the other dominated by activated pro- 
cesses. A first problem with the occurrence of these ad- 
ditional processes is that it hinders a clear separation 
of time scales between the fast transient dynamics and 
the structural relaxation, with which MCT exclusively 
deals, and thus makes the domain where the asymptotic 
regime of MCT is expected more difficult to resolve. In 
addition, the study of the corrections to the asymptotic 
laws shows that these laws are valid only in. a very small 
domain around the ideal glass transitionBQ: it can thus 
be expected that the dynamics in this domain will be 
markedly affected by crossover effects. Finally, from a 
more pragmatic point of view, since the asymptotic laws 
of the theory apply to rather specific aspects of the dy- 
namics (for instance, in the a relaxation regime, since 
the shape of the relaxation itself is non-universal, the 
main quantitative asymptotic prediction deals with the 
variation of the relaxation time with temperature), only 
selected parts of the available experimental information 
are usually involved in most experimental tests of the 
MCT on molecular glassformers, a situation which can 
be considered quite unsatisfactory. For these reasons, it 
appears desirable to have a means to test the MCT on 
experimental data which is not limited to the use of the 
critical predictions of the theory and which is able to give 
an overall account of the dynamics. This can precisely 
be achieved within the so-called schematic approach to 
experimental data. 

Schematic mode-coupling models are simple sets of 
mode-coupling equations retaining the essential non- 
linear structure of the equations of the MCT in one or 
two equations and with only a few parameters. They 
played a major role in the development of the theory to 
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investigate the universal properties of the general mode- 
coupling equation&.-aad— to, demonstrate the potentiali- 
ties of the theorylal3'tJ'yj'll3. Recently, they have been 
found useful as quantitative models for numerical fits 
of experimental dnta.fOicarmg a wide time or frequency 
rangeMyBBBBNM. The basic idea behind this 
schematic approach to experimental data is that, since 
all equations of the mode-coupling type share distinctive 
universal properties, meaningful informations are likely 
to be obtained if a suitable minimal set of equations of 
this family is used as an effective mode-coupling model to 
describe experimental data. One interest of the method 
is then that the universal properties of the mode-coupling 
equations are used in a milder sense than in their usual 
direct applications: it is thus possible to numerically take 
into account corrections to the asymptotic laws that can 
be important if the system is too far from its putative 
ideal glass transition. In addition, one can tentatively 
describe the full dynamics, including the non-universal 
a relaxation and the fast transient dynamics (although 
through rather crude approximations), and thus possibly 
overcome the problem of insufficient separation of time 
scales. For these reasons, the schematic approach ap- 
pears as a method of choice to critically investigate the 
validity of the previous tests of the theory based on the 
asymptotic laws. 

Of course, the power of the method comes at a price: 
even a minimal schematic model (in a sense to be defined 
latter) will have numerous adjustable parameters. This 
makes the approach numerically extremely flexible, and 
one finds for instance that it is possible to give an ac- 
count of the experimental dynamics with an ideal MCT 
schematic model even at low temperatures where ideal 
MCT is known inadequate. Thus one should not over- 



estimate the physical meaning of the fitted parameters. 
However, as was found in previous works, by critically 
looking at the results of a schematic analysis, one can ad- 
dress from a renewed point of view the problem of testing 
the predictions of the MCT on experimental data. 

In the present paper, the third of a series devoted to 
the schematic approach to experimental dataMltj, we aim 
to address various issues that have been overlooked in 
previous works, in order to better settle and understand 
the principles and results of this approach. In the next 
section, we give an overview of the schematic studies pro- 
posed so far and of the models on which they were based. 
In Sec. Ill, a set of simple requirements that need to be 
fulfilled in order to define a suitable minimal schematic 
model is proposed. 



[V, we examine the influ- 



In Sec. 

ence of the choice of the schematic model on the results 
of the analysis of a given set of experimental data, and 
Sec. ^ is devoted to a tentative validation of the method, 
by analyzing model spectra with well defined properties 
within the schematic approach. We finally present our 
conclusions and draw perspectives for the approach. 
II. MODELS FOR SCHEMATIC STUDIES 



In this section, we give an overview of the schematic 
studies proposed so far. This gives us the opportunity to 
settle a few definitions and notations to be used latter. 

The first schematic study has been proposed by Al- 
ba-Simionesco and KrauzmaixJ, where they analyzed 
low-frequency Raman scattering spectra of the fragile 
glassformer metatoluidine with various two-correlator 
schematic .models. One of them, based on the so-called 
F12 modelcJ'L3, is defined by the equations 



4>i(t) + rA&(t) + n?&(t) + o? / rn.it - T )4n(T)dT = o, % = o, i 

Jo 

m Q (t) = «i0o(f) + v 2 <j>o(t) 2 , 
m\(t) = rmo(t), 



(la) 

(lb) 
(1c) 



with initial conditions </>i(0) = 1, 4>i(0) = 0, i = 
0, 1, and has subsequently been used in studies of Ra- 
man spectra of other glassformersE-3 and, more impor- 
tantly, of depolarized light-scattering spectra of salol 
and CKN [CEI0.4K0. 6^03)1.4] covering a wide frequency 
range (more than four decades)lla. In all these works, 
where other choices for the jaemory-function mo have 
occasionally been consideredli3, the fitting function was 
defined, up to an amplitude factor, as the susceptibility 
associated with the effective correlator 

0schW = iMtf + (1 - 7)0i W 2 - (2) 
Another widely used model differs only from the pre- 



ceding one through the choice 

m l {t)=v s 4> {t)4> l {t), (3) 

an expression first proposed by Sjogren!^, and the use 
as the fitting function of the susceptibility associated 
to 0i, with an adjustable amplitude factor. It has 
been used to -study depolarized light-scattering spec- 
tra of glycerol!^ and orthoterphenyle (OTP)E3, as well 
as coherent neutron-scattering and Brillouin spectra of 
Nao.sLio.sPOsB An extended MCT version of this 
model has also been used recently to study depolarized 
light-scattering, dielectric-loss and incoherent neutron- 
scattering spectra of propylene carbonateEj. 
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Finally, self-intermediate scattering functions obtained 
by molecular dynamics simulations of ajaickcl-zirconium 
model have been analyzed by Teichlenia with a single- 
correlator model with mo defined as 



mo(t) = 1 (e^ o{t) 



1 



(4) 



At variance with the works mentioned above, a retarded 
damping term involving four adjustable parameters has 
been used in the present case instead of the simple Marko- 
vian approximation made in Eq. ( la ) involving the single 
parameter Tq. 

Looking at the equations, one sees that these mod- 
els involve a priori a large number of fitting parame- 
ters. To reduce this number, simplifying assumptions 
are usually made: in Refs. [l3[ [l4|, [l5|, and |l6|, the only 
temperature-dependent parameters are chosen to be v\ 
and V2, whereas in Ref. |l^, these are vi and v s , and in 
Ref. Uf, vi, v 2 and v s change with temperature, the latter 
depending also on the wave vector modulus. 



III. DEFINITION OF A SUITABLE MINIMAL 
SCHEMATIC MODEL 

In all schematic studies proposed so far, a certain de- 
gree of arbitrariness has often been involved in the choice 
of the details of the model to be used. The most serious 
case relates perhaps to the use of the memory function 
(Ik]) in our preceding worksOolij, which, at variance 
with Sjogren's modeled, cannot be given any precise the- 
oretical foundations. But other minor choices can also be 
questioned. There is indeed no physical reason to prefer 
using a simple correlator as the fitting function instead of 
a composite one like in Eq. (^|), the apparent advantage 
of the former choice, which seems to make the relation 
between the fitting function and the dynamical variable 
probed by the experiment transparent, relying on a mere 
tautology. One can also wonder why, for similar experi- 
mental data and schematic equations, the effective vertex 
entering the second memory function (r or v s ) has been 
held constant with temperature in Refs. [l^ and [l^ and 
not in Ref. 17 for instance, where it is v\ which has been 
kept fixed. In most cases, these various choices have been 
only validated a posteriori by the success of the model 
to fit the data under study and/or the obtention of pa- 
rameters whose variati ons .match the expectations from 
the microscopic theoryt§c3. 

There are nevertheless a few difficulties that could be 
expected, related to the possible existence of non-generic 
features of the chosen model. Indeed, as stated in the in- 
troduction, generic mode-coupling equations display re- 
markable universal propertiesu, but, when a schematic 
model is considered, because it is very simple, one can ex- 
perience severe limitations in the realization of these uni- 
versal features: certain values of important parameters of 
the theory can be unattainable or ad hoc correlations be- 
tween these parameters can be imposed. An example of 



the first situation is provided by the F n models, n > 2, 
for which 1 /2 is the only possible value of the exponent 
parameter A, whereas it can take any value in the inter- 
val [1/2, 1] for the general mode-coupling equations. The 
second situation is realized with the F12 model, where 
the exponent parameter and the critical non-ergodicity 
parameter f§ are always related through f§ = 1 — A. 

For these reasons, it seems desirable to understand 
what are the minimal requirements that must be met 
by a schematic model in order to avoid such biases and 
possible consequent inconsistencies in the results of its 
application. We propose here a series of criteria based 
on a discussion of the asymptotic properties of the mode- 
coupling equations in the vicinity of the ideal glass tran- 
sition. For this, we consider the case of an experiment 
probing the evolution of a normalized correlation func- 
tion <j)A associated to an unspecified dynamical process 
in a given glassforming liquid, when a thermodynamic pa- 
rameter X (temperature or density for instance) is varied. 
Note that the following discussion can easily be extended 
to situations where more tha«..©-ijie observable are consid- 
ered in a unified frameworloI3u3, with exactly the same 
conclusions. 

By looking at the asymptotic predictions of the the- 
orytl, a few important parameters emerge. On the one 
hand, one finds parameters that characterize the liquid 
under study and the thermodynamic path along which 
the experiment is performed. These are X c , the value of 
X at the ideal glass transition, Cx, the proportionality 
constant relating the separation parameter a defined by 
the theory to X through the relation 
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(5) 



and A, the exponent parameter governing numerous as- 
pects of the critical dynamics. On the other hand, one 
finds parameters that are specific to the correlation func- 
tion cj)A' f%, the non-ergodicity parameter at the transi- 
tion, and fiA, the so-called critical amplitude. Thus we 
identify five parameters that are expected to be of spe- 
cial relevance. The first three of them are involved in the 
description of the relative evolution of the two character- 
istic time scales of the dynamics in the ergodic phase, t a 
and tv characterizing the /3 regime and the a relaxation 
respectively. This can be seen from the ratio 



T -^ = B- 1 l h \Cx\- xl[2h) 



x — x r 



x,. 



-1/(26) 



(6) 



where b and B are uniquely determined by A§. In addi- 
tion, A alone determines the shape of the dynamics in the 
f3 regime. The two remaining parameters, j c A and Yia, are 
related to the amplitudes of the a relaxation and of the 
j3 process, as shown for instance by the scaling relation 



4>A{t) = f C A + h A \Cx\ 



1/2 



X - x r 



X r 



1/2 



g\{t/t a 



(7) 
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describing the (3 process and the onset of a relaxation in 
the ergodic phase close to the singularity. 

But in fact, this parameter counting is only valid when 
the full mode-coupling equations describing the dynamics 
of the system are known. Two differences appear when 
the only available information is the function <f>A at differ- 
ent X values. First, it can be seen that Ha never appears 
alone in the asymptotic laws, but always inside products 
h A \C x \ 1/2 - For this reason, the relevant parameter is 

1 /2 

rather the physical amplitude Ha = h A |Cx/(l — A)| ' , 
which appears for instance in the relation between X and 
the non-ergodicity parameter in the ideal glassy state, 



Sa = fl + H A 



X - X r 



x c 



1/2 



(8) 



This remark would be irrelevant if it was possible to know 
Cx from the data at hand, but this is not the case. In- 
deed, according to the theory, the shape of the a relax- 
ation is non-universal. This means that there exists an 
undetermined proportionality constant relating r a , the 
value emerging from the theoretical calculation, to the 
phenomenological a relaxation time r a (X), which can 
be defined for instance after 4>a{Jo) — V e and which 
is the only type of value that can be accessed from the 
knowledge of 4>a alone. Such a problem does not occur 
for the [3 process which is universal, i.e. it depends only 
on A. This means that the ratio (^) can only be known up 
to an undetermined factor and thus that the value of Cx 
cannot be obtained. This is in fact an important result, 
since it shows that one can know the absolute distance (in 
terms of the separation parameter a) of a system to its 
ideal glass transition only if the mode-coupling equations 
describing the full dynamics are known. 

Moving now back to the schematic approach to exper- 
imental data, one finds thus that, in order to analyze the 
evolution of the correlation function 4>a with an effective 
schematic correlator 4> sc h, the model must be able to re- 
produce the values of A, f A and Ha characterizing the 
system and the observable under investigation. A point 
on the transition surface of the schematic model must 
exist where one will have 
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(9) 



The last equality is not a strong constraint since h sc h and 
Cx,sch relate to decoupled aspects of the description of 
the dynamics, the first one being defined at the transition 
point, whereas the second one describes the approach of 
this point. One can thus expect that the interplay of 
these two aspects will allow to compensate possible non- 
generic features relative to /i sc h (like in the case of the 
F12 model where /i sc h = A sc h). The first two constraints 
are on the contrary very strong and need to be taken 
into account. One has in particular to worry about the 
accessible domain for the values of A sc h and / s c ch and to 



avoid the introduction of arbitrary correlations between 
these parameters. 

Eventually, one finds thus that a schematic model 
needs to have at least three adjustable coupling con- 
stants, at least one of them changing with temperature. 
Indeed, with a single coupling constant, only one pair 
(A S ch,/ s c c h) is possible. With two coupling constants, a 
relation linking A sc h and / s c ch exists, and for a given value 
of A sc h, only a finite number of / s c ch values will exist. One 
needs three coupling parameters to be able to decouple 
Asch and / s c ch . 

We can now look at the schematic models introduced 
in the previous section at the light of these new results. 
It is easy to see that these conditions are matched by the 
model used in Ref. [l^ and defined by Eqs. (|l|) and (j|), 
since it relies on four coupling parameters (v± , vi , r and 
7) to obtain the values of A sc h, / s c ch and J£. c h- This is 
not the case of the model used by TeichlercJ, which has 
only two coupling parameters 7 and n, and one indeed 
finds inconsistencies in the results: for instance, the expo- 
nent parameter A appears to change with the considered 
observable. It seems thus likely that the success of the 
calculation is essentially to be credited to the flexibility 
of the fitting procedure (in particular through the choice 
of a four-parameter retarded friction) and to the con- 
sideration of observables with high plateau values which 
probably turn out to be compatible with the value of 
the real exponent parameter. About the other studies, 
one finds in fact that the present discussion is inoperant, 
because only part of the dynamics is considered: the a 
relaxation is not present in Refs. |l3| and |l4|, and the fast 
transient dynamics is not included in Refs. |lq, [l7[ and 
with the effect of releasing in both cases the constraint 
on the value of / s c ch . 

An interesting consequence of these constraints is that 
they restrain the possibility to give a precise physi- 
cal meaning to the separate components of a schematic 
model. This is especially clear when a composite corre- 
lator like Eq. @ is used. Indeed, the relations between 
the components of the fitting function are expected to be 
mainly imposed by the need to reproduce the values of A, 
f A and Ha- it is then hardly conceivable that the value 
of the coupling constant 7 or the effect of the details of 
the expression for sc h will reflect any physical reality, 
for instance in terms of the light scattering mechanism. 
This is of course less obvious when the fitting function 
is a pure correlator, because of the tautological relation 
between the fitting function and the dynamical variable 
probed by the experiment, but, for the previous reasons, 
we believe that in general a schematic model has to be 
considered in its globality as a minimal phenomenologi- 
cal model allowing to give an account of the dynamics of 
a system in terms of mode-coupling effects, and nothing 
more has to be expected from it. 

To conclude this part, it is worth stressing that the 
present discussion refers only to the asymptotic predic- 
tions of the theory and that problems originating in pre- 
asymptotic corrections or the existence of activated pro- 
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cesses have been ignored. Moreover, only a few numbers 
have been considered. In particular, the precise shape of 
the relaxation (particularly in the a relaxation domain) 
has not been taken into account, whereas it can make 
an a priori well-suited model totally inefficient for the 
schematic approach. For these reasons, the preceding 
remarks are only guidelines aimed at rationalizing and 
optimizing the use of schematic models, but they offer 
no guarantee of success. 



IV. SCHEMATIC ANALYSIS OF 
EXPERIMENTAL SPECTRA WITH DIFFERENT 
MODELS 

In the present section, we propose a study of depolar- 
ized light-scattering spectra of a fragile glassformer by 
means of a variety of schematic models. The aim of this 
work is twofold. First, in all previous schematic analysis, 
only one model has ever been used at a time for a given 
set of experiments and no systematic investigation of the 
influence of the choice of the model on the results has 
been proposed so far. Secondly, in the previous section, 
a few parameters have been isolated that are expected to 
play a special role in the description of the dynamics of a 
given observable. The use of different schematic models 
to analyze the same set of experimental data, depending 
of their ability to provide consistent estimates for these 
values, should thus allow us to confirm or infirm the cor- 
responding discussion. 

The data we have studied are the CKN depolarized 
light-scattering spectra measured by Cummins et alx3 
and already analyzed with a schematic model in Ref. [l^. 
The reason CKN has been chosen (and not salol for in- 
stance) for this study is the low intensity of its a peak, 
which is slightly smaller than the microscopic peak: this 
guarantees that, when the spectra are studied with a 
composite effective correlator like in Eq. (0), the am- 
plitudes of the contributions of both terms will not be 
too different and thus this allows one to obtain a bet- 
ter contrast between the outcome of different schematic 
analysis. To limit the influence of the avoidance of the 
ideal glass transition, only rather high temperatures, be- 
tween 383 and 468 K, have been considered. Compared 
to our previous work, the two lower temperatures (363 
and 373 K) have been excluded, leading to small differ- 
ences in the results (within the error bars of the method) . 
Note that at the temperatures we have considered, one 
can expect the spectra to be free of the influence of the 
spurious experimental artifacts evidenced in Ref. |2^. 

Only schematic models closely related to the one used 
in Ref. |l5| and with the same number of parameters have 
been considered. As in this previous work, they have 
been used to perform simultaneous fits of the experimen- 
tal data at the nine selected temperatures: the whole dy- 
namical scenario developing as temperature varies is thus 
encoded in the only temperature dependent parameters 
of the schematic models, which are the effective vertices 
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FIG. 1: Experimental depolarized light-scattering spectra of 
CKN at 393 and 433 K (for readability, the data are reported 
every 4 points) and corresponding fitting curves obtained from 
a global fit of the nine spectra measured between 383 and 468 
K with two schematic models of type I, based respectively on 
the F12 and F13 memory functions. 



entering in the expression of the memory function mo- 

We first stick to the two-correlator F12 model. There 
are two routes by which the model can be modified. On 
the one hand, different second memory functions can be 
chosen. We have considered two cases: 

mi(t) = r(«i0o(*) + v 2 Mtf), (10a) 

mi(i)=«,0oW0i(*)- ( 10b ) 

On the other hand, one can vary the expression of the 
effective correlator as a function of </>o and <j>\. We have 
considered four possibilities: 



Wt) = 70o(*) 2 + (1-7)01 W 2 , 
f>sob(t) = 7 cAo(i) 2 + (1-7)01 (t), 
koh(t) = 70o(*) + (1-7)01 {tf, 
f>sob(t) = 70o (*) + (! -7)01 (*)■ 



(Ha) 
(lib) 
(11c) 
(lid) 



In previous worksEjEEl, the choice of a specific expression 
was originally motivated by assumptions on the light- 
scattering mechanism. However, from the discussion of 
the preceding section, it seems very unlikely that this 
question is of any relevance: indeed, since individually 
each correlator has no real physical meaning, the precise 
way it enters the calculation is clearly unimportant. 

A total of six models (numbered from I to VI, see Table 
H for definition) obtained by various combinations of the 
previous expressions have been used. Model I, defined by 
Eqs. (10a) and (|lla| ), was used in Ref. [If], whereas Model 
VI, using Eqs. (|I0b|) and (lid), is close to the one used in 
Refs .Jl6| , 17, and |18|. The fits are always very good (see 
Fig. Fyfor spectra and fits at two representative sample 
temperatures) and of equal quality for all models. The 
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FIG. 2: Effective vertices obtained from a schematic analysis of the CKN data with various F12 models. Left: temperature 
evolution of the effective vertices (open symbols: vi, closed symbols: Vi)- Right: trajectories in the control parameter plane 
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(lie) (He) (lid) 



r 

v s 
7 



8.6 3.6 5.9 2.2 



0.30 0.29 0.31 0.22 



14.2 
0.32 



7.5 
0.49 



T c (K) 388.5 388.9 388.5 
A 

t> ■ 

J sch 
^sch 
— Ct 
H sc h 



0.712 0.704 0.700 

0.448 0.452 0.446 

0.779 0.696 0.887 

0.150 0.183 0.113 

0.563 0.548 0.544 



388.8 
0.698 
0.449 
0.799 
0.133 
0.530 



390.1 
0.705 
0.424 
0.657 
0.230 
0.580 



390.8 
0.701 
0.429 
0.707 
0.150 
0.501 



TABLE I: Temperature-independent coupling constants and 
characteristic mode-coupling parameters obtained from a 
schematic analysis of the CKN data with various F12 models. 



corresponding results for the various coupling constants 
and relevant mode-coupling parameters are reported in 
Table | and Fig. |. 

Each model provides its own series of effective vertices 
v\ and V2 and values of the coupling constants r or v s 
and 7, distinct from those obtained with the other mod- 
els. The description of the dynamics in terms of effective 
mode-coupling parameters is thus clearly non-unique, as 
expected from the non-universality of the a relaxation. 
For all models, a behavior similar to the one found in 
Ref. [l5| is recovered, interpreted as follows: at the high- 
est temperatures the effective vertices are linear functions 



of temperature and show the progressive evolution of the 
system towards an ideal glass transition point, but, be- 
fore this point is reached, this behavior breaks down, 
evidencing the phenomenon of avoidance of the transi- 
tion and the inadequateness of the present ideal MCT 
approach at low temperature. The reader is referred to 
previous woxk for a discussion of the consequences of this 
observationtj. 

Looking at the trajectories in the (v±,V2) plane, one 
sees that the change in behavior of the effective vertices 
between the high and low temper ature domains is less 
visible when Sjogren's model, Eq. (10b), is used: in this 
case, the trajectories seem to behave rather regularly in 
the full temperature domain and to bend gently to fol- 
low the transition line, at variance with what is found 
from the use of Eq. (10a), where one can easily sepa- 
rate a domain where v\ does not vary much from a do- 
main where it decreases rather rapidly. We attribute this 
qualitative difference between the two models to the way 
the temperature dependence is entered in the schematic 
model at the level of the second memory function ra\. 
Indeed, at a given point (1)1,1)2) of the ergodic domain, 
one can qualitatively quantify the strength of the mode- 
coupling term in Eq. (|la|), i = 1, by considering the value 
of toi for representative values fo eauL.fi of the correla- 
tion functions <fio and <p\ , respectively!!!! When Eq. ( |l0b| ) 
is used with constant v s , the corresponding factor v s fof± 
depends only implicitly on temperature through the vari- 
ation of fo and fi with the position of the representa- 
tive point of the system in the control parameter space 
(vi,V2), whereas, when Eq. (10a) is used with constant 
r, the factor r(vifo + v 2fo 2 ) shows an additional explicit 
temperature dependence, since mi is defined directly in 
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terms of the coupling constants v\ and v%, and in fact 
this explicit dependence dominates the implicit one. If 
one now compares the evolution of these factors for both 
models when one explores the (i>i,i>2) plane, one eas- 
ily finds that, in order to have similar increases in the 
strength of the mode-coupling terms, one has to con- 
sider trajectories in the case of Sjogren's mod el th at lie 
above the corresponding trajectories for Eq. ( |lOa|) and 
that come rather close to the transition line, since in its 
vicinity only are found significant increases in / and /i, 
hence the difference in the results obtained with these two 
models: the catastrophic behavior of the vertices due to 
the avoidance of the transition is less visible in the case 
of Sjogren's model, because it is by construction already 
constrained to follow rather closely the transition line. 

Once this subtle non-generic feature has been recog- 
nized, one is able to rationalize the similarities and differ- 
ences between the results of the various schematic anal- 
ysis proposed so far. Since the light-scattering spectra 
of glycerol have been analyzed with the same parame- 
terization of Sjogren's model as the one used here for 
model V and VI, it is n£L surprise that the trajectory 
found by Franosch et alS3 is similar to those obtained 
here with these two models. This is the case of the 
study of Nao.5Lig.5PO3 as well, despite the choice of a 
temperature-dependent v s : above T c , the numerical fits 
result in a slow increase of this parameter when tempera- 
ture is decreased, which converts, as one would expect, to 
a trajectory in the (vi, 1)2) plane that is not different from 
what is obtained when v s is taken constant£3. The, case 
of OTP provides another illustration of this pointO. In- 
deed, in this work, the relative position of the vertex tra- 
jectory and the transition line is imposed from the start 
through the choice of a temperature-independent v\, a 
constraint that is compatible with the results obtained 
above for models I to IV, and, as could be anticipated 
from the previous discussion, this coincides with a rapidly 
increasing v s when temperature is decreased, i.e. a strong 
explicit temperature-dependence of the second memory 
function. The same seems to hold for the extended mode- 
coupling analysis of Gotze and VoigtmannE3, where a 
neat intersection of the ideal glass transition line by the 
trajectory in the parameter space of the schematic model 
is associated with values of v s increasing rapidly with de- 
creasing temperature. 

These observations appear at first sight rather anec- 
dotic. They have nevertheless at least one important 
physical implication. Indeed, comparing the results for 
glycerol and Nao.5Lio.5PO3, two liquids with intermedi- 
ate fragility, to those for salol, CKN, and OTP, that are 
all very fragile, it is tempting to propose a connection 
between the fragility of a glassformer and the shape of 
the effective vertex trajectory that describes its dynam- 
ics. The previous discussion rules out such a tentative 
connection and shows that the analogies and differences 
found between the various fluids simply originate in the 
details of the implementation of the schematic approach. 

We now turn to the quantitative analysis of our re- 



sults. For each schematic model, the values of the var- 
ious parameters characterizing the ideal glass transition 
have been obtained by extrapolation of the linear de- 
pendence of the effective vertices with temperature at 
high temperature, which we have identified as the sig- 
nature of the domain where the ideal MCT is adequate 
(in the next section, support will be given to this anal- 
ysis that differs from the one used in Refs. [l6] and fl8| ). 
They are reported in Table |. All six models provide 
consistent estimates for T c , A, / s c ch and H sc h- This con- 
sistency tends to validate the discussion proposed in the 
previous section about the relevant variables to describe 
a given observable in terms of mode-coupling parame- 
ters. This is clearly a non-trivial observation, since the 
values of the different parameters entering the calculation 
vary notably from one model to another (for instance, Ct 
changes by a factor of two between model III and V). 

In order to investigate stronger changes to the sche- 
matic model, we have considered other memory functions 
mo, keeping m\(t) = rm$(t) as in (10a). One of the roost 
studied schematic models after F12 is F13 defined asEj 



m (t) = vi<fi (t) + w 3 0o(i) 3 



(121 



We have studied two models based on this memory func- 
tion. They will be denoted I and IV, since they are di- 
rectly analogous to the preceding I and IV F i2 models. 
We observe a clear deterioration of the quality of the 
fits compared to F12, but it remains nevertheless accept- 
able (see Fig. [[]). In fact, the use of the F13 model is 
made difficult by one of its non-generic properties. In- 
deed, when a B-transition line has a termination point 
or when two transition lines meet, strong corrections to 
the generic asymptotic laws appear and noticeably affect 
the shape of the dynamicsE3E3E3. These corrections can 
make it impossible to use a schematic model to fit exper- 
imental data if the description of the system under study 
requires to come close to such singularities. In the F12 
model, a termination point is found at V\ = v-i = 1, but 
here A = I, which is quite far from the values character- 
izing the typical glassformers, of the order of 0.7. Thus 
no problem has to be expected a priori. But in the F13 
model, one finds that the A- and B-transition lines meet 
when i?i = 1, Vs = 4, with A = 3/4. This is rather close 
to the values usually associated to experimental systems 
and one can thus anticipate difficulties with this model, 
which we actually see. 

The results we have obtained with the F13 models are 
reported in Table [n| and in Fig. ||. The behavior of the 
effective vertices is the same as with F12, but slightly less 
regular. Critical parameters have been obtained from 
high temperatures as before. The agreement with the 
previous results for the values of T c and / s c ch is very good, 
but less good for the values of A and if S ch- However 
one can hardly conclude from these results, since these 
parameters characterize the f3 process, in which the non- 
generic limitations of the F13 model precisely show up. 

As it was stressed in the previous section, when using 
the schematic approach to experimental data, one has 
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model I IV 
r 1.79 0.62 
7 0.11 0.050 
T c (K) 389.3 389.7 
A 0.688 0.670 
/ s c ch 0.443 0.439 
Kch 0.541 0.444 
-C T 0.254 0.373 
H sc h 0.489 0.472 

TABLE II: Temperature-independent coupling constants 
and characteristic mode-coupling parameters obtained from 
a schematic analysis of the CKN data with two F13 models. 



to avoid any arbitrary correlation between important pa- 
rameters of the theory introduced in an ad hoc way by 
the choice of a specific schematic model. In all models 
studied so far, one can identify such a correlation in the 
fact that, along the B type transition line, an increase of 
/o is systematically associated to a decrease of A. This 
could have introduced a bias in the schematic calcula- 
tion that would have made the avoidance of the ideal 
glass transition line a purely numerical artifact of the 
method. Indeed, despite the absence of a true ergodicity 
breaking at T c in real systems, one expects that the in- 
crease of the non-ergodicity parameter below T c survives 
under the form of an increase of the a peak height of 
the susceptibility. The vertex drift found with the F12 
model could thus be simply traced back to the need to 
give an account of this evolution, the obtained trajec- 
tories corresponding precisely to an increase of f§ (and 



thus of / s c ch )ll3. It would then reflect a misconception and 
limitation of the schematic approach itself and not more 
generally a failure of the underlying ideal MCT as it was 
argued in Ref. 15. 

To show that our calculations are indeed free of this 
artifact, we have investigated a model in which the vari- 
ations of A and f§ can be decoupled. We chose the 
F29 model, the transition line of which is made of two 
branches along which A varies between 1/2 and 1: Be- 
tween the points v\ = 4, Vg = and u| = 27/7, 
Vg = 3 9 /(7 • 2 8 ), A and f§ grow simultaneously, whereas 
between the points v% = 27/7, w§ = 3 9 /(7-2 8 ) and v% = 0, 
Vg = 9 9 /8 8 , A decreases and f§ increases. Clearly, we 
have chosen this somewhat far-fetched model only for its 
mathematical properties and one could hardly give any 
physical motivation for the specific F29 functional form. 
However, we stress that, from the point of view of the 
criteria of the preceding section, this model is a priori 
free from any non-generic features and thus a perfectly 
acceptable schematic model: first, it does not show any 
limitation on the accessible values of A, at variance with 
all other F pq different from F12 i secondly, in the domain 
of A we are interested in (around 0.7), the influence of the 
higher-order A4 singularity from which its special math- 
ematical features originate is expected to be negligible. 

We have performed schematic analysis with this model 
in the vicinity of each branch of the transition line, in 
an attempt to find if the coupling of the evolutions of 
A and /q plays any role in our calculations. The cor- 
responding fits are quite poor: They only give a qual- 
itative account of the evolution of the dynamics and it 
would seem unreasonable to extract quantitative values 
for the critical parameters of the theory from the results. 
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However, a clear pattern manifests itself: Close to both 
transition branches, the behavior of the effective vertices 
corresponds to a decrease of the exponent parameter, ir- 
respective of the evolution of /q. In fact, this could have 
been anticipated from the fact that, below the value of 
T c we have determined, the maximum of the a peak lies 
outside the available frequency window, and this elimi- 
nates the requirement to reproduce its amplitude at these 
low temperatures. However it is reassuring to show by 
an independent calculation that our results are free from 
this potential bias, that should certainly be taken into ac- 
count if data where available on a wider frequency range. 

To summarize, we have shown in this section that dif- 
ferent schematic models applied to a same set of exper- 
imental data provide qualitatively similar results. We 
have found that they are in good quantitative agreement 
for the determination of the values of the parameters 
that have been distinguished in the previous section, pro- 
vided no non-generic features, the importance of which 
has been emphasized for the F13 model, interfere with 
their determination. 



V. SCHEMATIC ANALYSIS OF MODEL 
SPECTRA 

In the previous section, we have tried to validate the 
schematic approach and its results showing that differ- 
ent models provide consistent results when applied to 
a same set of experimental data. However, on various 
aspects, this study is not quite satisfactory. Indeed, de- 
spite favorable indications gained from the use of the F29 
model, the meaningfulness of the anomalous behavior of 
the effective vertices at low temperature remains to be 
assessed. In our previous worwla, it has been attributed 
to the existence of relaxation processes that are not taken 
into account in the ideal MCT and that make it impossi- 
ble to describe consistently the dynamics of the studied 
systems with a schematic model based on this theory. To 
strengthen this point of view, it would be necessary to 
show that the schematic approach is actually able to give 
an account of a complex purely mode-coupling dynam- 
ics with smoothly varying fitting parameters, despite of 
the extreme simplicity of the schematic equations com- 
pared to the complexity of the equations of the original 
problem. Moreover, it would be interesting to analyze 



data for which the values of the essential parameters A, 
f A et Ha are known in advance, in order to check to 
what extent the schematic approach is able to provide 
quantitative estimates for them. 

A possibility to perform such a test is to calculate se- 
ries of curves from mode-coupling equations known in 
advance and to analyze them as one would do for any 
set of experimental data. It remains then to compare 
the results obtained within the schematic approach with 
the properties that can be directly determined from the 
starting equations. We report the results of such a test 
in the present section. 

Recently, Franosch et al. have proposed a detailed 
study of the MCT predictions for the collective dynamics 
of a simple liquid, with special emphasis on the asymp- 
totic laws valid in the vicinity of the ideal glass transi- 
tion and their correctionsO. Because their model has been 
very well and extensively characterized, it appears clearly 
as the model of choice to be the basis of the present eval- 
uation of the schematic approach. Besides it has already 
been used in a similar spirit for a critical study of the 
method consisting of the use of numerical fits of the (3 
correlator for experimental tests of the MCTtEl. 

The model is based on the mode-coupling equations for 
the Brownian dynamics of the hard-sphere fluid with par- 
ticle diameter d, the thermodynamic state of which is de- 
termined by a single parameter, the compacity 77, related 
to the number density p through the relation rj = irpd 3 /6. 
One starts with the generalized Langevin equations for 
the time evolution of the Fourier components of the nor- 
malized density-density correlation functions (f> g , 

T q <]) q + (f)q+ / drm q (t - r)0 3 (r) = 0, (13a) 
Jo 

with the initial condition a ^ g (0) = 1. T q is chosen equal 
to t m i c S q /(qd) 1 , where t mlc is some arbitrary molecular 
time scale and S q is the structure factor of the fluid. The 
mode-coupling memory kernel, defined as a double inte- 
gral in Fourier space (thus a six-dimensional integral), 
is approximated after moving to bipolar coordinates, by 
the Riemann sum of the resulting two-dimensional inte- 
gral on a grid of M equidistant values of the wave vector 
modulus, with a grid step h such that the qd values vary 
between h/2 and (2M — l)h/2. One obtains 



ph 3 S M min ^ M ' 
m i( t ) = ^EzhJ2 J2 kf>S P S k 



32#7T 2 <f 



(P + q 2 - p 2 )c k + (f 



k 2 )c p 4>k{t)4> p (t), 



(13b) 



k=l p=|q-fc| + l 



with qd = qh, kd = kh, pd = ph and q, k, p £ tion function related to the structure factor through the 
{1/2,3/2,... , (2Af — l)/2}. c q is the direct correla- Ornstein-Zernike relation pc q = 1 — l/S q . As in pre- 
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vious worksS, the static quantities have been computed 
within the Percus-Yevick approximation, M = 100 and 
h = 0.4 have been used, and the length and time units 
have been chosen such that d = 1 and t m ; c = 160. One 
finds then that the ideal glass transition of the model 
occurs at a compacity r/ c comprised between 0.51591213 
and 0.51591214. 

This model aims to represent the dynamics of a hard- 
sphere colloidal suspension and has been shown to com- 
pare very favorably at a quantitative level with the ex- 
perimental data of van Megen et a£3. 

The equations (13) have been solved following the al- 
gorithm described in Rcf. 23 , then the associated suscep- 



tibilities have been computed by Fourier transform using 
Filon's algorithm. We concentrate on the susceptibilities, 
for which it has been shown that the range of validity of 
the asymptotic laws is usually smaller than for the cor- 
responding correlation functionstj, in order to be as close 
as possible to the situations in which the schematic ap- 
proach is usually applied to experimental data. 

Five values of the compacity such that the distances 
to the ideal glass transition are s n = (r/ c — i] n )/r/ c = 
10""/ 3 , n = 4,5,6,7,8 have been considered. Larger 
n values have been excluded, since they would corre- 
spond to spectra covering a much wider frequency range 
compared to the experimental situation. We have ana- 
lyzed spectra for three values of the wave vector modulus: 
qi = 7.0, corresponding to the main peak of the struc- 
ture factor, q2 — 10.2, corresponding to the first mini- 
mum after the main peak, and = 14.6, where certain 
pre-asymptotic corrections seem to be rather small. 

Having now our reference spectra to be used as input 
like experimental data would be, each series of curves at a 
given q value has been fitted with the two correlator F12 
model referred to as model IV in the previous section. 
For compatibility reasons, no inertial term is included, 
i.e. one has 

Ti<]>i + <pi+ [ drrmit - t)^(t) = 0, i = 0,l.(14) 
Jo 

Note that an inertial term would be required both in the 
microscopic theory and in a schematic model to study 
a molecular glassformer, the dynamics of which is New- 
tonian, and to reproduce its oscillatory dynamics in the 
Raman band or in the boson-peak region. However, for 
the goal we wish to achieve here, i.e. an assessment of the 
shape of the vertex trajectory and the calculation of the 
critical parameters of the theory from model theoretical 
spectra, the fact that it has been consistently omitted 
in the original model and in the schematic equations is 
irrelevant. No amplitude factor is needed, since the data 
are already normalized (to check this point, a few calcu- 
lations have been done with such a factor, showing that 
it remains equal to 1). As above, only the effective ver- 
tices V\ and Vi are allowed to vary with <j). This means 
in particular that constant values of tq and t± are im- 
posed, despite the fact that the r q values of the original 
equations actually change with compacity through S q . 
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FIG. 4: Susceptibility spectra of a model simple liquid (sym- 
bols) and fits obtained with an F12 schematic model (contin- 
uous lines), for n = 8 and three values of the wave vector. 
The size of the symbols has been chosen to roughly model the 
typical dispersion of data points obtained from experiments. 
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0.102 




0.329 




0.503 




C 


0.5153 


0.5159 


0.5155 


0.5159 


0.5154 


0.5159 


A 


0.678 


0.735 


0.676 


0.735 


0.679 


0.735 


rc 
Jq 


0.846 


0.849 


0.414 


0.419 


0.268 


0.269 


hq 


0.306 


0.323 


0.738 


0.642 


0.603 


0.582 


Co 


2.48 


1.54 


1.36 


1.54 


1.81 


1.54 


H q 


0.849 


0.779 


1.513 


1.548 


1.433 


1.403 



TABLE III: Comparison of the critical parameter values ob- 
tained from a schematic analysis of model spectra and from 
a direct study of Eqs. (jL3|) . 



The quality of the fits is good, comparable to what 
is obtained when experimental light-scattering spectra 
are analyzed, as illustrated in Fig. |^. Values obtained 
for the various adjustable mode-coupling parameters are 
reported in Table III and in Fig. 

As usual, the effective vertices are the important pa- 
rameters of the calculation. To each series of curves cor- 
responds a different set of effective vertices, stressing the 
fact that the study of a single observable is not enough 
to give an unequivocal account of the dynamics of the 
system. Here again, this is a consequence of the non- 
universality of the a relaxation as described in the MCT 
framework. The important observation is that the behav- 
ior of the effective vertices is perfectly regular: v\ and 
are approximately linear with e (and thus v± is linear 
with V2). Nothing similar to what has been found with 
experimental data close to the transition line shows up 
in the results. This tends to confirm the assumption that 
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the anomalous behavior found from experiments and at- 
tributed to the avoidance of the ideal glass transition is 
not connected to a technical limitation of the schematic 
approach, but actually, at a more fundamental level, to 
the inadequateness of ideal MCT in a temperature do- 
main starting noticeably above the extrapolated ideal 
glass transition point. 

We can now determine the critical parameter values 
by extrapolating the linear behavior we have found and 
compare the results with those obtained directly from 
the original equations, in order to see if the procedure we 
have used previously is meaningful and if the correspond- 
ing results are quantitative. The values are reported in 
Table [II. Of the three parameters whose importance 



has been stressed in Section |III| is clearly the one 



for which the agreement between the original value and 
the schematic calculation is the best. The determina- 
tion of H q is less precise but remains satisfactory. On 
the contrary, a notable disagreement on the value of the 
crucial exponent parameter A is found. The difference 
we find is significant and shows unambiguously that the 
schematic approach does not allow a priori to obtain the 
exact value of the exponent parameter, at least when the 
studied dynamics only covers a frequency range of the 
order of four or five decades. From this point of view, 
the schematic approach suffers from the same difficulties 
as the use of the (3 correlator, which has been shown nnt 
to be able to give the value of A in the same conditionsEl 
It has however to be stressed that, even if it is not per- 
fect, the schematic approach turns out to perform much 
better than the latter approach. Surprisingly, it is found 
that the three series of curves do not provide values scat- 
tered around the exact value A = 0.735, but instead, the 



same underestimated value is consistently found at all 
three wave numbers. It seems thus very likely that no 
improvement would be obtained if the three observables 
were analyzed in a unified framewo*kp4n the spirit of the 
most recent schematic calculationsll§E3. 

The origin of the present problem is strictly connected 
to the existence in the analyzed data of corrections to 
the asymptotic behavior predicted by the theory and not 
to a deficiency of the schematic approach. Indeed, when 
spectra corresponding to larger values of n (i.e. corre- 
sponding to compacity values closer to the critical one) 
are analyzed, a slight change in the behavior of the ef- 
fective vertices occurs around n — 10 and the trajectory 
described in the (v\, 112) plane bends to indicate the point 
of the transition line corresponding to the correct value of 
A. We stress, to avoid any ambiguity, that this change of 
behavior has nothing to do with the one found from ex- 
perimental data: its amplitude is very small and it occurs 
in a domain where the points representative of the system 
in the (v\, V2) plane are extremely close to the ideal glass 
transition line, meaning that the dynamics spreads over 
more than eight decades in frequency and corresponding 
for a real glassformer like CKN to a regime where the be- 
havior of the effective vertices is already strongly affected 
by the existence of activated processes. 



VI. CONCLUSION AND PERSPECTIVES 

In this Paper we have proposed a detailed study of the 
schematic mode-coupling approach to experimental data. 
Our aim was twofold. 

First, we wanted to settle the whole approach on firmer 
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theoretical grounds, through a better understanding of 
the parameterization it involves and a clarification of the 
nature of the results it can provide. Indeed, the schematic 
approach involves a certain degree of arbitrariness in the 
choice of the models to be used and the introduction of 
numerous adjustable parameters, which are never con- 
sidered as good things when it comes to testing a the- 
ory. Thus in order to credibilize the whole approach, 
one has to address these a priori problematic issues and 
show that they are actually harmless. Hence this in- 
vestigation, whose outcome is the formulation of simple 
rules that need to be taken into account in order to de- 
fine minimal schematic models suitable for use in this 
context and the isolation of a few parameters which are 
particularly relevant for the description of the dynam- 
ics and for which the schematic approach should pro- 
vide estimates. Following these c riteria ,, it rttpys out that 
most previous schematic studiesHoEjtj'EjH are based 
on well designed schematic models and use parameteri- 
zations which are close to the minimal one, in particular 
by letting only effective vertices vary with temperature. 

Secondly, our goal was a validation of the approach, 
through the investigation of its robustness and of its po- 
tential quantitative character by a variety of numerical 
tests. The results are rather satisfactory. Indeed, we 
have been able to validate the results obtained in Ref. [l5| 
and the way they were analyzed. We have shown that 
they have a limited, if any, dependency in the schematic 
model we use and we have been able to identify and rule 
out many potential biases that could have affected our 
results. Most properties of the dynamics whose rele- 
vance have been emphasized in the theoretical discussion 
of the approach appear to be quantitatively evaluated, 
with the notable exception of the exponent parameter A, 
for which the method is not too bad, but still far from 
perfect. It has thus to be stressed that the precise and 



unambiguous determination of this parameter, which ac- 
tually plays a very important role in the theory, remains 
elusive, by any of the methods that are available so far: 
all calculations indeed seem to be too sensitive to the 
presence of significant pre-asymptotic corrections (as is 
shown here or, in the case of the use of the j3 corre- 
lator, by Ref. [33]) and very probably to the existence of 
crossover effects due to-the non-ideal behavior of real and 
simulated glassformersEJ. 

In conclusion, we believe that the present work paves 
the way to new applications of the schematic approach. 
Indeed, the interest of the method to analyze experimen- 
tal data has been emphasized in all previous works and 
we think its relevance has been further strengthened here. 
However, at this level, it was only a kind of simple phe- 
nomenological mode-coupling approach and the main is- 
sue was its ability to consistently describe a set of ex- 
perimental data with a simple model. Now, from our 
results, it seems possible to go one step further and to 
use this approach as an effective means to obtain precise 
estimates of transition temperatures, non-ergodicity pa- 
rameters and critical amplitudes to be used for instance 
in quantitative tests of the MCT on simple systems that 
are amenable to a full first-principle treatment. 
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